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ABSTRACT 

We compare the CO (J =1-0) and HI emission in the Large Magellanic Cloud (LMC) in 
three dimensions, i.e. including a velocity axis in addition to the two spatial axes, with the 
aim of elucidating the physical connection between giant molecular clouds (GMCs) and their 
surrounding HI gas. The CO J =1-0 dataset is from the second NANTEN CO survey and the 
HI dataset is from the merged Australia Telescope Compact Array (ATCA) and Parkes Telescope 
surveys. The major findings of our analysis are: 1) GMCs are associated with an envelope of 
HI emission, 2) in GMCs [average CO intensity] oc [average HI intensity] 1 - 1 * ' 1 and 3) the HI 
intensity tends to increase with the star formation activity within GMCs, from Type I to Type 
III. An analysis of the HI envelopes associated with GMCs shows that their average linewidth is 
14 km s _1 and the mean density in the envelope is 10 cm -3 . We argue that the HI envelopes are 
gravitationally bound by GMCs. These findings are consistent with a continual increase in the 
mass of GMCs via HI accretion at an accretion rate of 0.05 Moyr -1 over a time scale of 10 Myr. 
The growth of GMCs is terminated via dissipative ionization and/or stellar- wind disruption in 
the final stage of GMC evolution. 

Subject headings: Magellanic Clouds — ISM: clouds — ISM: molecules — ISM: atoms — galaxies: 
ISM — radio lines: ISM 



1. Introduction 

Giant molecular clouds (GMCs), the most massive aggregations of interstellar matter with 10 5_6 M©, 
are the principal sites of star formation in galaxies. It is important to understand how GMCs are formed 
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out of the less dense atomic interstellar gas in order to understand galactic evolution. The interstellar HI 
gas has densities of less than several 10 cm~ 3 while molecular clouds have densities larger than 100 cm~ 3 . 
It is reasonable to assume that HI is being converted into H 2 either by thermal/gravitational instabilities 
and/or shock compressions, although the detailed processes of this conversion arc not yet well-understood. 
Sato & Fukui (1978) and Hasegawa, Sato, & Fukui (1983) identified cold HI gas associated with GMCs in 
M17 and W3/4 and suggested that the cold HI gas may be converted into molecular gas for these GMCs. 
Subsequently, Wannier et al (1983) showed that five molecular clouds are associated with warm HI envelopes 
and suggested that such HI envelopes may be common around GMCs. Nonetheless, associations between 
GMCs and HI envelopes are difficult to identify systematically throughout the Galactic disk, since GMCs 
samples are restricted to the solar vicinity due to the crowding effects (Andersson, Wannier, & Morris 1991). 
As a consequence, the CMC-HI association has not been well established. 

The Magellanic system - including the LMC, the SMC and the Bridge - is an ideal laboratory to study 
star formation and molecular cloud evolution because of its proximity to the Milky Way (e.g., Fukui et al 
1999; Mizuno et al. 2001; Mizuno et al. 2006; Fukui et al. 2008; Ott et al. 2008; Kawamura et al, 2009). 
We expect that the Magellanic system can also shed light on the physical connection between GMCs and 
their atomic surroundings. Indeed, the LMC may offer the best place for such a study because of its nearly 
face-on orientation and level of star formation activity. The LMC's molecular cloud population, which is 
best traced via the CO emission, provides a key to understand the galaxy's star formation. Molecular clouds 
are able to highlight the location of star formation due to their highly clumped distribution in both space 
and velocity. The LMC's atomic gas, by contrast, has lower densities and is only weakly coupled to sites of 
active star formation, but it is the most promising candidate for the mass reservoir of GMC formation (e.g., 
Blitz et al. 2007). We note, moreover, that cold HI gas has been detected in the LMC (e.g., Dickey et al. 
1994), and that the correlation between CO and HI may provide crucial observational evidence about the 
molecular cloud formation process. 

Wong et al. (2009) compared the HI and CO emission throughout the LMC on a pixel-by-pixel basis 
using the second NANTEN CO and ATCA+Parkes HI datasets. These authors studied correlations between 
the integrated CO and HI intensities, where the latter was integrated over all velocities with HI emission 
or over individual Gaussian components. They found that CO emission is associated with high intensity HI 
gas but that intense HI emission is not always associated with CO. They also discovered a weak tendency 
for CO to be associated with HI components that have relatively low velocity dispersion. This suggests that 
energy dissipation of the HI gas may be required for the formation of molecular clouds. Following the global 
analysis by Wong et al. (2009), we focus here on the HI associated with individual GMCs in the LMC. 
In order to address this issue, we conduct a detailed comparison between the CO and HI emission in three 
dimensions, i.e. (x,y, v), at a spatial resolution of ^40 pc and a velocity resolution of 1.7 km ss _1 . The 
present study is complementary to the work by Wong et al. (2009); in conjunction, the two studies provide 
a new insight into the CO-HI connection. In section 2, we briefly review the basic observational properties 
of GMCs in the LMC. In section 3, we describe our method of analysing the CO-HI correlation and present 
our results. We discuss the physical interpretation of our results in section 4 and provide a summary of our 
major conclusions in section 5. 



2. GMCs in the LMC; the second NANTEN CO survey 



The LMC is extended by more than 30 square degrees across the sky and it has been a difficult task to 
make a systematic survey of the CO emission at angular resolutions sufficient to resolve individual GMCs. 
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Fukui et al. (1999) made such a survey in the 2.6mm CO line with the NANTEN 4m mm-wave telescope 
and published the first results in Fukui et al. (1999). Subsequently, these authors completed another survey 
of the LMC, improving the sensitivity by a factor of two. The second NANTEN CO survey has cataloged 
272 GMCs (Fukui et al. 2008). The basic physical parameters of GMCs in the LMC are similar to those 
in the Milky Way and other nearby galaxies. Their masses range from ~ 10 5 M Q to ~ 10 6 M Q ; the mass 
spectrum is quite steep with a slope of dN/dM ~ M~ 2 . The X factor, the ratio of the H 2 column density 
to CO intensity, is - 7 x 10 20 cmT 2 (K km s" 1 ) ~ x (Fukui ct al. 2008; Blitz et al. 2007). The complete 
sampling of the NANTEN survey has also allowed us to make a statistical study of GMCs with various young 
objects including HII regions and young stellar clusters. Kawamura et al. (2009) confirmed that there are 
three classes of GMCs, that can be categorised according to their association with young star clusters as 
originally indicated in the analysis of the first NANTEN LMC survey (Fukui et al. 1999; Yamaguchi et al. 
2001). Type I GMCs show no signs of active star formation in the sense that no O stars are being formed. 
Type II GMCs are associated with small HII region(s), indicating the formation of isolated O stars, but do 
not host any stellar clusters identified by Bica et al. (1996). Type III GMCs are actively forming stars as 
shown by their association with large HII regions and young stellar clusters. These classes are interpreted 
as an evolutionary sequence from Type I to III; the lifetime of a GMC is estimated to be a few x 10 Myr 
in total (Kawamura et al. 2009). The stage after Type III is likely to be a very violent dissipation of the 
GMC due to the UV photons and stellar winds produced by the nascent clusters, as seen spectacularly in 
the region of 30 Dor (Yamaguchi et al. 2001). A comparison of the physical parameters of the GMCs shows 
that the size and mass of the clouds tend to increase from Type I/II to Type III. A summary of the average 
mass and size for the three GMC classes is presented in Table 1, as taken from Kawamura et al. (2009). 



3. Correlation between CO and HI 
3.1. 3-D correlation 

Previous studies of star formation in galaxies have employed 2-dimcnsional (2-D) maps of HI intensity 
with large spatial averaging on scales between ~100 pc and 1 kpc (e.g., Schmidt 1972; Kennicut et al. 1988). 
Here, we make a 3-D comparison between the CO and HI in the LMC where the 3-D datacubes have a 
velocity axis in addition to two spatial axes projected on the sky. Preliminary results of the comparison 
have been published elsewhere (Fukui 2007). We use the 3-D datacube of CO obtained with NANTEN 
(Fukui et al. 2008) and an HI datacube obtained with ATCA and Parkcs (Kim et al. 2003). The CO 
emission traces GMCs and the HI emission traces less dense atomic gas. Figure 1 shows an overlay of the 
velocity-integrated intensities of the CO and HI emission; from this, it is clear that GMCs in the LMC tend 
to be located towards HI filaments or local HI peaks, suggesting that HI is a prerequisite for GMC formation 
(Blitz et al. 2007). However, it is also clear that there are many HI peaks and filaments without CO emission 
(Wong et al. 2009). Figure 2 shows typical CO and HI line profiles in the LMC. The CO emission is highly 
localized in velocity: the HI emission ranges over 100 km s _1 while the CO emission has a typical linewidth 
of less than 10 km s^ 1 . We note that the large velocity dispersion of HI may be dominated by physically 
unrelated velocity components along the line of sight, i.e. the HI gas associated with the GMC may only 
be the small fraction of HI with velocities close to that of the CO emission. Previous studies of the CO-HI 
connection that use velocity integrated 2-D maps may therefore overestimate the intensity of the associated 
HI emission along each line of sight. By making use of the velocity dimension, the present 3-D analysis may 
allow us to identify the HI gas that is physically connected to the GMCs. The NANTEN and ATCA+Parkes 
datacubes have somewhat different spatial and velocity resolutions, so we have convolved both datasets to 
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a spatial resolution of 40 pc x 40 pc, and a velocity resolution of 1.7 km s . The total number of 3-D 
pixels is approximately 2 x 10 6 across the area of surveyed by NANTEN. The HI and CO intensities are 
expressed in units of ?b(K) and (K); the 3er noise levels of the HI and CO datacubes are 7.2 K and 
0.21 K respectively. Strictly speaking, the one-to-one correspondence between a velocity and a position is 
not guaranteed because there is a chance that physically unrelated HI gas may have the same velocity as 
HI gas related to the GMC along the same line of sight. Our results identify the HI associated with GMCs 
and suggest that such contamination along individual sightlines may not be a serious problem. We further 
note that HI absorption towards background radio continuum sources does not affect significantly the HI 
intensity at the present spatial resolution, as verified toward 30 Dor, one of the brightest radio continuum 
sources in the LMC. 

Figure 3 shows a histogram of the HI intensity in the 3-D datacube. Pixels with the significant CO 
emission (T^ > 0.21 K) are shown in red. Throughout this paper, histograms always use the values of 
3-D pixels with no spatial integration unless otherwise stated. The histogram in Figure 3 shows that the 
fraction of CO-detected pixels increases monotonically with the HI intensity, suggesting the HI intensity is a 
necessary condition to form GMCs, consistent with the conclusion by Wong et al. (2009). About one third 
of the pixels with Tb(HI) of ~ 90 K exhibit CO emission, but it seems that there is no sharp threshold value 
of HI intensity that is required for GMC formation. 

Figure 4 shows a histogram of the HI intensity for the three GMC types. Each pixel detected in CO 
belongs to one of the GMCs cataloged in Fukui et al (2008). Figure 4 clearly shows that the HI intensity 
tends to increase from Type I to Type III, although the dispersion is considerable. The average HI intensity 
for Types I, II and III GMCs across the LMC is 34 ± 16 K(lcr), 47 ± 17 K and 56 ± 19 K, i.e. the average HI 
intensity increases with the level of star formation activity within the GMC. In order to test for variation 
within the galaxy, we tentatively divide the galaxy into three regions, i.e., Bar, North and Arc, as shown in 
Figure 4(right). Histograms for each region, shown in the lower three panels of Figure 4, reveal the same 
trend, suggesting that the present trend is common over the whole LMC. The total number of the pixels in 
each region is 429,510, 666,060, and 405,093. Type I, Type II and Type III GMCs include 330, 1,065, and 
639 pixels in the Bar region; 346, 878, and 1,158 pixels in the North region; and 389, 650, and 1,231 pixels 
in the Arc region. 

Figure 5 shows the velocity channel maps of the HI distribution associated with Type I, Type II, and 
Type III GMCs. These channel maps clearly show that the HI is associated with the CO. The CO distribution 
has small structures of ~ 100 pc or less and the HI appears to be associated with the GMC on larger scales of 
~ 100 to 400 pc. The HI emission is not always symmetric with respect to a GMC, even though HI typically 
envelopes each GMC. The associated HI is often elongated along the GMCs and the region of intense HI 
emission is usually < 100 pc wide. The CO emission typically extends over a velocity range of ~ 5 km 
s _1 ; beyond a few times this velocity range, the associated HI emission generally becomes much weaker or 
disappears. 



3.2. Physical properties of the HI envelope 

In general, it is a complicated task to derive reliable physical properties of the HI gas associated with a 
GMC because the HI profiles are a blend of several different components along the line of sight, making it 
difficult to select the HI gas that is physically connected to a GMC. Another obstacle is that the HI emission 
is spatially more extended than the CO emission and has a less clear boundary than the CO. 
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For our analysis, wc first selected GMCs with simple single-peaked HI profiles from the Fukui et al 
(2008) catalog. The resulting sample consists of 123 GMCs in total. Their catalog numbers and basic 
physical properties, taken from Fukui et al. (2008), are listed in Table 2. For these GMCs, we tested 
whether there was a bias in their location with respect to the kinematic center of the galaxy, in their CO 
lincwidth or in their molecular mass. The histograms in Figure 6 indicate that there is no particular trend for 
these properties of the selected GMCs compared to GMCs in the complete catalog, suggesting that there is 
no appreciable selection bias. We applied a Kolmogorov-Smirnov test to the three histograms and calculated 
maximum deviations of 0.031, 0.061 and 0.117 respectively for the three parameters. These values are less 
than the critical deviation, 0.129, for a conventional significance level of 0.05, confirming that there is no 
selection bias. 

Next, we made Gaussian fits to the HI and CO profiles towards the CO peak of each GMC. This 
procedure yields a peak intensity, peak velocity and half-power linewidth for each line profile (a summary 
is given for each GMC type in Table 1). Figure 7 shows the relation between the CO lincwidth and the 
difference between the CO and HI peak velocities. We find the HI and CO peak velocities to be in good 
agreement, showing only a small scatter of less than a few km s" 1 . Figure 8 shows two histograms of the 
HI and CO linewidths. We see that the HI linewidth is typically 14 km s _1 , roughly three times larger 
than that of CO. Figure 9 shows a correlation between HI and CO linewidths. The two quantities show a 
positive correlation with a correlation coefficient of 0.39. The correlation coefficient is determined using the 
Spearman rank method throughout the present paper. The kinematic properties of HI and CO, as illustrated 
in Figures 7 and 9, lend further support to a physical association between the HI and CO. 

In order to estimate the size of the HI envelope surrounding each GMC, we construct an HI integrated 
intensity map of each GMC. First, we find the local peak in the HI intensity cube surrounding the CO 
emission, and then integrate the HI intensity over the velocity channels corresponding to the FWHM of the 
HI line profile at this peak position. Next we estimate the area, S, where the HI integrated intensity is greater 
than 80 % of the value at the local HI peak. We then calculate the radius of the HI envelope, i?(HI) from 
its projected area, S = 7ri?(HI) 2 . The HI integrated intensity is calculated for all the pixels with detectable 
CO emission; the spatial distribution of the HI emission generally shows a peak and a reasonably defined 
boundary. The 80 % level was chosen after a few trials using different levels; it is the maximum value for 
which a reasonable HI size is obtained for 116 of the 123 envelopes. While 80 % seems to be rather high for 
such a definition of a cloud envelope, the HI size can be unrealistically large compared to the CO cloud size 
along a filamentary HI distribution if we use a lower level such as 60 %; an example is shown in Figure 5b. 
The HI radius is then corrected for beam dilution by adopting Gaussian de-convolution with an HI FWHM 
of 2.6 arcmin, i.e. the same procedure that is applied to the radius of the GMC determined from the CO 
emission, R(CO) (Kawamura et al 2009). R(CO), the CO peak position, i?(HI), the HI peak position, and 
the deviation between the peaks in parsecs, \J (Act 2 + AS 2 ), are listed in Table 2. Seven clouds, for which 
R(HI) comprises only a few pixels, are denoted by asterisks. The physical extent of the HI envelopes can be 
as large as a few 100 pc. Ws long as we use the local HI peak of individual GMCs as our reference point, it 
seems difficult to contrive an alternative uniform definition of the envelope size is difficult using the present 
HI dataset. 

In Figure 10, we show a histogram of the spatial deviation of the HI and CO peaks. This shows that 
nearly 80 % of the HI envelopes peak within 120 pc of the local CO peak, and that nearly 60 % of the HI 
envelopes peak within 80 pc of the CO peak. These separations may seem large compared to i?(HI), but we 
argue that they are reasonable if the HI is enveloping CO at scales over <~ 300 pc. It should be noted that 
the HI envelopes are not concentric with the CO emission, but are rather "enveloping" with some offsets in 
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peak positions as illustrated in Figure 5. We thus expect to find some difference in general between the peak 
positions of the CO and HI emission (as seen in Figure 10), but the fact that the majority of HI peaks are 
located within 120 pc of the CO peaks is clearly suggestive of a physical association between the GMCs and 
their surrounding atomic gas. In Figure 11, we show a correlation between i?(HI) and R(CO) for 62 GMCs 
whose radius is greater than 30 pc and find that they are positively correlated with a correlation coefficient 
of 0.45. Despite the relatively flatter distribution of the HI, the size of the HI envelope does seem to correlate 
with the size of the GMC. 

To summarize our analysis in this section, we find for the 123 GMCs with single-peaked HI profiles that: 
1) the peak velocities of the CO and HI are in good agreement (Figure 7); 2) the CO and HI linewidths 
show a positive correlation (Figure 9); 3) the HI envelopes, defined using the 80 % level of the local HI 
integrated intensity peak, are mostly (~ 80 %) centred within 120 pc of the peak CO position (Figure 10); 
and 4) the radius of the HI envelope is positively correlated with the size of the GMC for GMCs with radii 
greater than 30 pc (Figure 11). These four results lend further support to the idea that HI envelopes are 
physically associated with GMCs, reinforcing the impression conveyed by a global comparison between the 
HI CO emission in the LMC (Figure 3) and the morphological similarity between the CO and HI in individual 
velocity channels. (Figure 5). 

Next, we made an estimate of the HI column density for the 123 GMCs by using the relation (N (HI) 
[cm -2 ] = 1.8 x 10 18 / T h dv[K km s -1 ]). The average values for the three GMC types are listed in Table 
1. We find the peak HI column density is mostly in a range of 2-5 xlO 21 cm~ 2 . We estimate the typical 
density in the HI envelopes to be ~ 10 cm -1 by dividing the peak HI column density 2-5 xlO 21 cm~ 2 by the 
typical size of the associated HI 50-100 pc (see Figure 11). The mass of the HI envelopes is large, typically 
~ 5 x 10 4 M Q . The HI envelopes are likely gravitationally bound by GMCs because a half of the HI line 
width, 7 km s _1 , is nearly equal to y/GM/R ~ 6 km s _1 for M = 2 x 10 5 M Q and R = 40 pc, the average 
values of Type II GMCs. 

In Figure 12, we plot the relationship between the average CO and HI luminosity of the 123 GMCs. 
For each GMC, we selected pixels where CO emission is significantly detected: only these pixels are used 
to calculate the HI and CO luminosity of each cloud. In order to derive the average CO luminosity of a 
GMC, we estimated /(CO) K km s _1 by dividing the sum of all the significant CO emission within the GMC 
by the projected area of the GMC. /(HI) K km s _1 is calculated in a similar manner, summing up the HI 
emission from the pixels within the GMC with significant CO emission and dividing the result by the area 
of the GMC. The regression shown in Figure 12 is well fitted by a power law with an index of 1.1, indicating 
a nearly linear correlation between /(CO) and /(HI) in a GMC. 

4. Discussion 

4.1. GMCs with HI envelopes 

The present analysis has successfully identified the HI envelopes associated with GMCs on the basis of 
a 3-D analysis of GMCs in the LMC. The HI intensity in the envelope depends on the star- forming activity 
within the GMC in the sense that the integrated HI intensity in the envelope increases from Type I to Type 
III (Figure 4, section 3.1). In other words, massive GMCs have massive HI envelopes and less massive GMCs 
have less massive HI envelopes. 

The HI intensity is a product of the spin temperature and the optical depth of the HI 21 cm transition, 



-7- 



provided that the line is optically thin. This is likely to be the case as the HI profiles show few hints of 
saturation like a flat top. The observed maximum HI brightness temperature is around 100 K and this 
suggests that the HI spin temperature is significantly higher than 100 K. Therefore we infer that the HI 
intensity should represent optical depth and, accordingly, HI column density, if the spin temperature is 
roughly uniform across the LMC. The HI spin temperature in the LMC may be higher than in the Galaxy 
due to a more intense UV field and the lower dust extinction (Israel 1996; Imara & Blitz 2007; Dobashi et 
al. 2008). For the sake of discussion, we shall tentatively assume that the spin temperature lies between 
150-600 K and is fairly uniform in the HI envelope. We note that the HI mass is accurately determined 
as long as the HI emission is optically thin. The HI mass does not depend on Ts under the optically thin 
assumption, because the level populations of the spin doublet having only ~ 10~ 5 eV is well thermalized in 
any realistic density range due to the slow magnetic dipole decay in ~ 10 7 yrs. 

Fukui et al. (1999) suggest that the three classes of GMC indicate an evolutionary sequence from Type 
I to Type III in a few 10 Myr (instead of "Type", these authors used "Class" with the same meaning). 
Kawamura et al. (2009) present a more detailed analysis of the association between GMCs and young stellar 
clusters, confirming Fukui et al.'s evolutionary scheme. These studies indicate that only the youngest star 
clusters with an age less than ~ 10 Myr are clearly associated with GMCs, and that older clusters with 
an age greater than 10 Myr are not associated with GMCs. Assuming a steady state scenario, this implies 
that the natal gas of clusters is quickly disrupted within 10 Myr. Considering the complete sampling of 
both clusters and GMCs in the LMC, this strongly suggests that the population of Type III GMCs must be 
replenished on timescales of 10 Myr. Since the typical time scale of GMC formation is at least 10 Myr, as 
estimated by the crossing timescale - i.e., the cloud size divided by its velocity dispersion, 100 pc/10 km/s 
= 10 Myr, a measure of the minimum timescale for GMC formation - we expect to have a similar population 
of Type III GMCs and Type III precursors. A straightforward interpretation is that Type I and Type II 
GMCs are these precursors (see for details Kawamura et al. 2009). An alternative possibility is a more 
ad-hoc situation in which Type III GMCs are formed suddenly in a few Myr by an external disturbance, 
such as a dynamical interaction. Such a strongly time-dependent scenario seems unlikely, however, since the 
three GMC types are fairly uniformly distributed over the LMC (Kawamura et al. 2009). Figure 4 of this 
paper also shows that the three classes of GMCs are distributed across the galaxy. 

We have also seen that Type III GMC tend to be more massive than Type I and Type II GMCs (Table 2; 
Kawamura et al. 2009). A natural interpretation within the evolutionary scenario is that the enveloping HI 
gas accretes onto a GMC and that the GMC mass increases with time. The accreted HI envelope is converted 
into H 2 in ~ 10 Myr due to increased density, optical extinction and UV shielding. This infall scenario is 
consistent with the linear relationship between /(HI) and /(CO) in Figure 12; by contrast, growth of GMCs 
via collisions between HI clouds would have a steeper relationship, with /(CO) proportional to /(HI) 2 . For 
an infall scenario, the infall motion can arise from the gravity of a GMC and possibly from a converging flow 
driven by super bubbles, while the thermal motion is negligibly small 1.4-3 km s _1 for kinetic temperatures 
of <~ 150-600 K. We can roughly estimate the infall velocity to be half of the HI linewidth, i.e. ~ 7 km s _1 . 
This value is consistent with the free fall velocity, ~ 6 km s _1 , for a typical Type II GMC. For spherical 
accretion, where the GMC is surrounded by an HI envelope with a radius of ^40 pc, volume density of n(HI) 
^10 cm~ 3 , and an infall speed of ~ 7 km s _1 , we estimate the mass accretion rate to be ~ 0.05 M Q yr _1 . 
Over the typical timescale of the GMC evolution, i.e. <~ 10 Myr, the increase in molecular mass amounts to 
~5x 10 5 M Q , which is roughly consistent with the observed typical value for the mass of a Type III GMC 
(~4x 10 5 M Q , Table 1). In the evolutionary picture, the mass accretion of a GMC is terminated by the 
violent disruption and/or ionization of the molecular material by stellar winds and ionization from young 
stars. 
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The infall scenario offers a reasonable interpretation of the HI and CO properties of GMCs that we 
have explored in this paper. It remains to be seen, however, if an infall velocity field is consistent with 
2-D observations; careful analysis of an isolated HI envelope with little kinematical disturbance or nearby 
contamination could be used to verify this. It is also important to clarify whether the HI linewidth is affected 
by turbulence to a significant degree. 

It could be argued that the HI gas surrounding GMCs is supplied by the recombination of HII into 
HI, as both Type II and Type III GMCs are associated with HII regions. This alternative seems unlikely, 
however: first, the HII regions in Type II GMCs are compact and therefore do not constitute a significant 
mass reservoir; second, the HI envelope in Type III GMCs are not spatially well-matched with the HII 
regions and clusters. In N159, for instance, the HII regions and young clusters are confined to the north of 
the GMC, wherease the HI is more widely distributed in the east and south (Figure 13). 

4.2. HI-H 2 conversion in GMCs 

This study has shown that the HI and CO distributions correlate well on 40-100 pc scales. It is worth 
noting, however, that the correlation becomes less clear on smaller scales of ~ 10 pc within a GMC. Figure 
13 shows an overlay of the HI and CO distributions for the Type III GMC N159, where the CO data was 
obtained using the ASTE 10 m sub- mm telescope in the 12 CO J =3-2 emission line (Mizuno et al., 2008). 
Figure 13 shows that the HI becomes less bright at T b (HI) 70-80 K toward N159E at (R.A., Dec.)=(5 h 40 m , 
—69°, 45'), compared to intensities of Tb(HI) ~ 120 K in the HI envelope. A similar behavior was noted by 
Ott et al. (2008). This is unlikely to be cuased by absorption of the radio continuum emission, as there is 
no radio continuum emission toward N 159E. We regard this behaviour to be illustrative of the conversion 
of HI into H 2 , as well as the the lower spin temperature in the interior of a GMC. In the inner part of a 
GMC, HI is converted into H 2 via reactions on grain surfaces on a time scale of ~10 Myr. The HI density 
is typically ~1 cm~ 3 , compared to a total molecular density of a few 100 cm~ 3 , corresponding to atomic to 
molecular hydrogen ratio of -100 (e.g., Allen & Robinson 1977; Spitzer 1978; Goldsmith et al. 2007). The 
spin temperature is also lower, and is likely to equal the molecular gas kinetic temperature of — 60 K (e.g., 
Sato & Fukui 1978; Mizuno et al. 2009). In the HI envelope, on the other hand, the spin temperature is 
probably between — 150-600 K and the HI density is estimated to be ~ 10 cm -3 with no H 2 . The lower 
lb (HI) in the interior is likely due to the lower spin temperature and the lower HI density. The mass of the 
apparently cold HI gas towards N159E is approximately 10 % of the mass of the HI envelope, ~ 10 5 M Q , 
if we assume that the cold HI is optically thin, which suggests that the cold HI within the GMC is not a 
dominant mass component of the atomic+molecular cloud complex. 

It has been shown that there are cold HI components in the LMC as measured from emission and 
absorption observations toward radio continuum sources (Dickey et al. 1994). These authors detected 
HI absorption features toward 19 of 30 continuum sources in the LMC and argued that T s of the cold 
components can be as low as 40 K. Such cold HI components may be associated with GMCs. It is however 
not clear observationally how the cold HI in absorption is related to GMCs because none of the absorption 
measurements by Dickey et al (1994) coincide with the NANTEN GMCs. 

An issue which has been raised in Wong et al. (2009) is that higher HI intensity is a necessary but not 
sufficient condition for CO formation. In other words, there are many places with high HI intensities in the 
LMC without CO. The current analysis has focussed solely on the HI gas surrounding GMCs and therefore 
does not directly address this issue. We note, however, that it is possible that HI gas with the same intensity 
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can have a significantly different density. We propose an interpretation that the atomic gas in regions with 
high HI intensities but no CO may have lower densities and higher temperatures. This interpretation could 
be tested in the future by high velocity-resolution HI observations that can resolve subtle variations in HI 
profiles and hence identify spatial variations of the atomic gas temperature. 

5. Summary 

We have carried out the first 3-D analysis of the connection between the CO and HI emission in a galaxy. 
The major results of our study are as follows: 

1. A 3-D comparison at a resolution of 40 pc x 40 pc x 1.7 km s™ 1 has revealed that the fraction of HI 
associated with CO tends to increase as a monotonic function of HI intensity without a sharp threshold 
for the CO formation. 

2. We find that GMCs are associated with HI envelopes on scales of ~ 50-100 pc. The HI envelopes have 
typical volume densities of ~ 10 cm -3 and an average linewidth of ~14 km s -1 , which is about three 
times larger than the linewidth of CO. We argue that the HI envelopes are gravitationally bound by 
GMCs. 

3. For 123 GMCs with single-peaked HI profiles, we find a correlation such that [average CO intensity] oc 
[average HI intensity] 11 * 01 . There is a clear increase of the associated HI intensity from GMC Type 
I to Type III. 

4. We interpret our results to mean that a GMC increases in mass via continuous HI accretion over 
a timescale of ~10 Myr and with a mass accretion rate of 0.05 M Q yr _1 , before being disrupted by 
ionization and stellar winds from young clusters. The accreted HI is likely to be converted to molecular 
hydrogen due to the higher shielding within a GMC. 
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by a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and 
Technology of Japan (No. 15071203), from JSPS (No. 14102003, No. 18684003, and core-to-core program 
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Tabic 1 . Physical properties of GMCs 



GMC Typo 


Number of 




i? a 


N(HI) b 


AVlsr (HI) b 


AVlsb. (CO) a 




GMCs 


(xlO 5 M ) 


(pc) 


(xlO 21 cm ~ 2 ) 


(km s" 1 ) 


(km s" 1 ) 


Typo I 


72 


2(2) 


37 (16) 






5.0 (2.5) 


Type II 


142 


2 (3) 


33 (19) 






4.8 (2.2) 


Type III 


58 


5 (10) 


51 (36) 






6.9 (3.0) 


Type I (selected) 


24 


2 (3) 


35 (17) 


2.4 (0.9) 


13.9 (4.0) 


4.5 (2.1) 


Type II (selected) 


67 


2 (3) 


41 (22) 


2.6 (1.2) 


14.6 (4.1) 


4.4 (1.6) 


Type III (selected) 


32 


4 (3) 


55 (23) 


3.3 (1.5) 


16.1 (3.3) 


5.5 (1.5) 



a Fukui ct al. (2008); Kawamura ct al. (2009) 

b Half intensity full width derived by gaussian fitting. 

c Sclcctcd clouds having single peaked HI profiles. 



Note. - 



Average propeties of the GMCs. The values in parentheses are the standard deviation. 



Tabic 2. List of selected 123 GMCs 









peak position(CO) c 


i?(CO) a 


peak position (HI) c 


_R(HI) d 


N(UI) C 


VAa 2 + A<5 2f 




numbcr a 


name 3. 


TvDe b 




o \ ij lyou j 


he] 




oyo lyou ) 


foci 


10 21 [cm- 2 ] 


foci 


comment^ 


i 


LMC N J0447-6910 


I 


4h47.7m 


-69° 14' 


44 


4h47.7m 


-69° 14' 


34 


2.6 







4 


LMC N J0449-6910 


III 


4h49.1m 


-69°16' 


72 


4h49.5m 


-69° 14' 


74 


2.7 


91 




5 


LMC N J0449-6826 


II 


4h49.5m 


-68°28' 


100 


4h50.1m 


-68°32' 


77 


1.7 


152 




9 


LMC N J0450-6930 


11 


4h50.5m 


-69°36' 


33 


4h50.2m 


-69°34' 


55 


2.1 


83 




11 


LMC N J0451-6858 


I 


4h51.3m 


-69°4' 


29 


4h51.3m 


-69°4' 


29 


2.7 







12 


LMC N J0451-6704 


III 


4h51.9m 


-67° 8' 


108 


4h52.2m 


-67°8' 


62 


3.2 


70 




17 


LMC N J0453-6909 


III 


4h54.ini 


-69°14' 


85 


4h53.3m 


-69° 12' 


59 


2.1 


159 




22 


LMC N J0455-6830 


II 


4h55.9m 


-68° 36' 


47 


4h55.9m 


-68°36' 


18 


2.2 







23 


LMC N J0455-6634 


III 


4h56.4m 


-66° 42' 


64 


4h56.4m 


-66°40' 


* 


1.9 


28 


30,32,33,236 


24 


LMC N J0455-6930 


III 


4h56.3m 


-69°36' 


29 


4h56.3m 


-69°36' 


18 


0.8 







26 


LMC N J0457-6844 


III 


4h57.5m 


-68°48' 


29 


4h57.5m 


-68°48' 


24 


1.5 







27 


LMC N J0457-6826 


III 


4h57.4m 


-68° 30' 


40 


4h57.ini 


-68°32' 


18 


1.1 


91 




35 


LMC N J0459-6614 


II 


4h58.8m 


-66° 18' 


79 


4h58.4m 


-66°24' 


85 


2.7 


130 


30,32,33,236 


36 


LMC N J0500-6622 


11 


5h0.7m 


-66°26' 


52 


5h0.7m 


-66° 26' 


37 


1.5 







38 


LMC N J0502-6903 


II 


5h2.4m 


-69°6' 


37 


5h2.4m 


-69°6' 


41 


2.1 







39 


LMC N J0503-6553 


11 


5h3m 


-65°58' 


40 


5h3m 


-65°58' 


62 


2.7 







40 


LMC N J0503-6540 


II 


5h3.1m 


-65°44' 


37 


5h3.4m 


-65°46' 


24 


0.6 


75 




41 


LMC N J0503-6828 


II 


5h3.6m 


-68° 32' 


33 


5h4.3m 


-68°36' 


37 


1.1 


160 




43 


LMC N J0503-6643 


11 


5h3.4m 


-66°48' 


33 


5h3.4m 


-66°50' 


24 


1.2 


26 




45 


LMC N J0503-6719 


III 


5h3.7m 


-67° 24' 


44 


5h3.7m 


-67° 24' 


24 


2.0 







46 


LMC N J0504-6802 


III 


5h4.6m 


-68°6' 


47 


5h4.3m 


-68°8' 


18 


1.0 


83 




48 


LMC N J0504-7007 


II 


5h5m 


-70° 12' 


44 


5h5.1m 


-70°12' 


8 


1.4 







49 


LMC N J0504-7056 


III 


5h5.3m 


-71°0' 


52 


5h5.3m 


-71°0' 


34 


1.0 







52 


LMC N J0506-6753 


II 


5h6.8m 


-67° 58' 


23 


5h7.2m 


-67° 56' 


24 


1.0 


83 




53 


LMC N J0507-7041 


II 


5h7.9m 


-70° 46' 


23 


5h7.8m 


-70°48' 


24 


0.9 


28 




55 


LMC N J0507-6858 


I 


5h8.1m 


-69°2' 


47 


5h8.1m 


-69°0' 


29 


1.3 


28 




57 


LMC N J0508-6905 


I 


5h9.1m 


-69°8' 


23 


5h9.2m 


-60°8' 


8 


1.0 







60 


LMC N J0509-6827 


II 


5h9.4m 


-68° 32' 


33 


5h9.5m 


-68°30' 


47 


1.6 


26 




61 


LMC N J0509-7049 


II 


5hl0.2m 


-70° 54' 


29 


5h9.8m 


-70° 54' 


8 


1.0 


87 




62 


LMC N J0509-6912 


I 


5hl0.2m 


-69° 16' 


52 


5hl0.6m 


-69°16' 


37 


1.7 


79 




63 


LMC N J0510-6853 


III 


5hl0.4m 


-68° 56' 


84 


5hl0m 


-68° 52' 


77 


2.3 


93 


245,68 


64 


LMC N J0510-6706 


II 


5hl0.7m 


-67° 12' 


52 


5hll.lm 


-66° 58' 


74 


1.2 


222 


67 


65 


LMC N J0511-6927 


11 


5hll.6m 


-69° 30' 


23 


5hl2m 


-69°32' 


47 


0.8 


83 


244,71 


67 


LMC N J0512-6710 


II 


5hl2.ini 


-67° 14' 


55 


5hl2.4m 


-67° 14' 


55 


1.6 


79 


64 


68 


LMC N J0512-6903 


II 


5hl2.2m 


-69°4' 


80 


5hll.4m 


-69°4' 


88 


1.8 


166 


63 



Table 2 — Continued 









jjcdiv puaii 


IU111 KjKJ J 


-R(CO) a 


pLdI\ JJUSllj 


iuii 1 1 ii i 


-R(HI) d 


7V(HI)° 


\l An 2 4- AiS 2f 


11 LI III U CI 


a 


TVnp b 
±ype 


n>( R1 Q^fA 

Lx. 1 _D i_y OU J 


o \-D j. you j 


LP C J 


tjtyn> j_ you j 


uyiD i_y ou j 


fnrl 

LP C J 


-1U [CII1 


IP C J 


69 


LMC N J0512-7028 


II 


5hl2.8m 


-70°32' 


68 


5hl2.8m 


-70°30' 


50 


1.5 


26 


71 


LMC N J0513-6936 


II 


5hl3.5m 


-69° 40' 


79 


5hl3.5m 


-69° 44' 


59 


1.9 


62 


72 


LMC N J0513-6922 


III 


5hl3.5m 


-69° 26' 


49 


5hl3.2m 


-69° 24' 


34 


1.7 


83 


74 


LMC N J0514-7010 


II 


5hl4.9m 


-70° 14' 


66 


5hl4.9m 


-70° 14' 


34 


1.6 





77 


LMC N J0515-7034 


I 


5hl5.6m 


-70° 38' 


44 


5hl5.6m 


-70°38' 


18 


0.6 





80 


LMC N J0516-6807 


II 


5hl6.1m 


-68° 2' 


124 


5hl7.2m 


-68°4' 


62 


1.4 


237 


81 


LMC N J0516-6616 


I 


5hl6m 


-66°20' 


29 


5hl6m 


-66° 18' 


24 


0.3 


26 


83 


LMC N J0516-6922 


11 


5hl7m 


-69° 26' 


23 


5hl7.4m 


-69°26' 


34 


1.3 


87 


84 


LMC N J0516-6559 


II 


5hl6.7m 


-66° 2' 


29 


5hl7.1m 


-66°2' 


29 


0.6 


70 


86 


LMC N J0517-7114 


III 


5hl8.3m 


-71° 18' 


44 


5hl8.3m 


-71°20' 


37 


1.3 


26 


89 


LMC N J0517-6642 


II 


5hl7.6m 


-66° 44' 


40 


5hl7.6m 


-66°44' 


18 


0.9 





90 


LMC N J0517-6932 


II 


5hl8.5m 


-69°36' 


33 


5hl8.5m 


-69°36' 


* 


0.9 





92 


LMC N J0518-7001 


I 


5hl8.8m 


-70°6' 


33 


5hl8.1m 


-70°4' 


18 


1.1 


168 


93 


LMC N J0518-6620 


I 


5hl8.3m 


-66° 24' 


23 


5hl8m 


-66°22' 


* 


0.5 


75 


94 


LMC N J0518-6951 


II 


5hl9.2m 


-69°54' 


23 


5hl9.2m 


-69°56' 


24 


0.9 


26 


95 


LMC N J0519-6625 


I 


5hl9m 


-66° 28' 


23 


5hl9m 


-66°30' 


* 


1.2 


26 


96 


LMC N J0519-6938 


III 


5hl9.6m 


-69°40' 


72 


5hl9.6m 


-69° 40' 


29 


2.0 





97 


LMC N J0519-7113 


I 


5h20m 


-71° 18' 


40 


5h20.4m 


-71°16' 


24 


0.9 


91 


99 


LMC N J0520-7043 


I 


5h20.4m 


-70° 44' 


52 


5h20m 


-70°46' 


37 


1.0 


91 


100 


LMC N J0520-665 


II 


5h20m 


-66° 54' 


37 


5h21m 


-66° 50' 


52 


1.3 


226 


103 


LMC N J0521-714 


II 


5h21.7m 


-71°46' 


23 


5h21.7m 


-71°46' 


24 


1.2 





104 


LMC N J0521-701 


II 


5h21.6m 


-70° 16' 


49 


5h21.2m 


-70° 18' 


29 


0.8 


91 


105 


LMC N J0521-700 


II 


5h21.6m 


-70°4' 


79 


5h21.6m 


-70°6' 


44 


1.4 


26 


106 


LMC N J0521-684 


II 


5h21.5m 


-68° 44' 


29 


5h21.5m 


-68°44' 


8 


1.7 





108 


LMC N J0521-684 


11 


5h21.8m 


-68°50' 


55 


5h21.8m 


-68°50' 


34 


1.4 





110 


LMC N J0522-694 


III 


5h22.7m 


-69° 44' 


52 


5h22.3m 


-69°44' 


18 


1.2 


79 


115 


LMC N J0522-654 


III 


5h22.6m 


-65° 44' 


57 


5h22.6m 


-65°42' 


44 


1.7 


26 


118 


LMC N J0523-682 


I 


5h22.9m 


-68° 24' 


44 


5h22.9m 


-68° 26' 


37 


3.2 


26 


119 


LMC N J0523-664 


II 


5h23m 


-66° 46' 


57 


5h22.7m 


-66° 46' 


41 


2.5 


79 


123 


LMC N J0523-713 


11 


5h24.2m 


-71°42' 


49 


5h26.4m 


-71°38' 


62 


1.4 


466 


126 


LMC N J0524-672 


II 


5h24.5m 


-67° 30' 


44 


5h24.5m 


-67° 30' 


24 


2.3 





127 


LMC N J0524-702 


I 


5h24.8m 


-70°30' 


37 


5h24.4m 


-70°30' 


24 


0.6 


87 


130 


LMC N J0524-691 


II 


5h24.9m 


-69° 18' 


33 


5h24.5m 


-69° 12' 


47 


1.2 


123 


131 


LMC N J0524-713 


11 


5h25.5m 


-71°36' 


29 


5h25.9m 


-71°38' 


64 


1.7 


99 


132 


LMC N J0525-694 


II 


5h26.2m 


-69°52' 


33 


5h25.8m 


-69°56' 


34 


1.0 


99 



comment 6 



91 



100 



131,249 



123,249 
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peak posit 


ion(CO) c 


i?(CO) a 


peak posit 


ion(HI) c 


i?(HI) d 


7V(HI)° 


VAa 2 + A<5 2f 




number a 




name* 


Type b 


a(B1950) 


<5(£1950) 


[PC] 


a(B1950) 


5(_B1950) 


[PC] 


10 21 [cm- 2 ] 


[PC] 


comment 6 


133 


LMC 


N 


J0525-O91 


III 


5h26m 


—69° 20' 


33 


5h24.5m 


—69° 12' 


47 


0.8 


359 


130 


134 


LMC 


IN 


J 0525-062 


II 


5h25.7m 


—66° 24' 


23 


5h25.3m 


— 66°26' 


116 


3.8 


75 


135 


136 


LMC 


IN 


J 0526-083 


II 


5h26.2m 


—68° 38' 


37 


5h26.2m 


— 68°38' 


41 


1.4 







139 


LMC 


M 
IN 


JUDZ0-0o4 


III 


5h26.6m 


—68° 50' 


29 


5h26.3m 


— 68°48' 


29 


1.9 


91 




141 


LMC 


IN 


J 0526-055 


II 


5h26.5m 


—65° 56' 


23 


5h26.2m 


— 66°0' 


160 


2.0 


93 


134,135,250 


143 


LMC 


IN 


J0526-711 


II 


5h27.5m 


— 71°22' 


97 


5h26m 


— 71°26' 


120 


1.3 


346 


153,156 


145 


LMC 


M 
IN 


TO ^07 7HQ 
JUOZ 1- l\JO 


II 


5h28m 


Tr* 000/ 

—70 38 


40 


5h28.1m 


—70 40 


18 


0.6 


26 




146 


LMC 


M 
IN 


JUOZ 1- /Uo 


II 


5h28.6m 


—70° 54' 


37 


5h28.6m 


—70 54 


72 


1.0 





149 


149 


LMC 


IN 


J 0528- (05 


II 


5h29m 


— 71°0' 


33 


5h29.4m 


— 71°0' 


24 


1.7 


87 




150 


LMC 


IN 


J 0529-083 


II 


5h30.2m 


— 68°32' 


29 


5h30.6m 


— 68°36' 


142 


1.8 


106 


154,163 


153 


LMC 


IN 


J0530- / 10 


III 


5h31.6m 


—71° 10' 


80 


5h31.6m 


—71° 10' 


47 


2.6 







154 


LMC 


M 
IN 


JUool-Ooo 


III 


5h31.7m 


—68 34 


108 


5h32.7m 


— 68°28' 


55 


2.7 


243 


150,163 


155 


LMC 


IN 


J0532-O/4 


III 


5h32m 


—67° 46' 


68 


5h32.7m 


—67° 46' 


44 


3.1 


157 




156 


LMC 


M 
IN 


JUooz- (11 


II 


5h32.9m 


—71° 16' 


47 


5h32.9m 


—71° 16' 


41 


2.5 







157 


LMC 


IN 


JUooz-Doo 


II 


5h32.5m 


—68 40 


40 


5h32.5m 


n a r\f 

—68 40 


41 


2.7 







158 


LMC 


IN 


J 0532-062 


III 


5h32.4m 


—66° 30' 


47 


5h32m 


— 66°26' 


29 


0.9 


99 




162 


LMC 


IN 


Tnc^Qo cot: 
J0532-O85 


II 


5h33.4m 


— 69°0' 


33 


5h33m 


— 69°0' 


18 


2.4 


79 




171 


LMC 


IN 


J 0535-090 


III 


5h36m 


— 69°4' 


66 


5h36.4m 


— 69°4' 


91 


3.5 


79 


183 


172 


LMC 


IN 


J 0535-084 


II 


5h36.2m 


—68 46 


37 


5h35.8m 


—68 46 


41 


1.6 


79 




179 


LMC 


M 
IN 


JUDO /-DDI 


III 


5h37.3m 


—66° 20' 


47 


5h37.3m 


— 66°20' 


44 


2.1 







180 


LMC 


IN 


J 053 (-062 


II 


5h37.4m 


—66° 28' 


23 


5h37.1m 


— 66°30' 


34 


1.8 


75 




184 


LMC 


N 


J0538-693 


II 


5h38.7m 


—69° 36' 


33 


5h39.1m 


— 69°34' 


41 


3.5 


83 


191 


190 


LMC 


N 


J0538-685 


III 


5h38.9m 


—68 56 


23 


5h41.1m 


—68° 54' 


81 


2.9 


480 




191 


LMC 


N 


J0539-693 


III 


5h39.4m 


— 69°32' 


57 


5h39.1m 


— 69°34' 


41 


3.5 


83 


184 


205 


LMC 


N 


J0542-711 


III 


5h42.ini 


-71°20' 


64 


5h42m 


-71°18' 


47 


2.5 


28 




207 


LMC 


N 


J0542-694 


II 


5h43.5m 


-69°46' 


23 


5h43.5m 


-69° 46' 


18 


4.1 





266 


208 


LMC 


N 


J0543-675 


III 


5h43.5m 


-67° 58' 


37 


5h43.8m 


-67° 58' 


34 


2.0 


79 


211 


209 


LMC 


N 


J0543-692 


II 


5h43.9m 


-69° 28' 


33 


5h44.2m 


-69° 22' 


47 


3.0 


106 


214,216,268 


211 


LMC 


N 


J0543-675 


III 


5h43.8m 


-67° 56' 


23 


5h43.8m 


-67° 58' 


29 


1.8 


28 


208 


214 


LMC 


N 


J0543-691 


II 


5h43.4m 


-69° 16' 


33 


5h44.2m 


-69°22' 


44 


2.5 


211 


209,216,268 


215 


LMC 


N 


J0544-712 


I 


5h45.2m 


-71°28' 


59 


5h44.7m 


-71°28' 


57 


1.9 


87 




218 


LMC 


N 


J0544-671 


II 


5h44.9m 


-67° 20' 


29 


5h44.4m 


-67° 28' 


50 


1.0 


167 


213 


220 


LMC 


N 


J0545-694 


III 


5h45.9m 


-69° 50' 


47 


5h45.6m 


-69°52' 


37 


2.0 


83 




222 


LMC 


N 


J0546-710 


I 


5h46.8m 


-71°6' 


44 


5h46.8m 


-71°6' 


44 


1.7 







223 


LMC 


N 


J0546-693 


II 


5h46.8m 


-69° 38' 


52 


5h47.2m 


-69°40' 


57 


3.3 


99 





Table 2 — Continued 









peak posit 


ion(CO) c 


.R(CO) a 


peak posit 


ion(HI) c 


i?(HI) d 


iV(HI) e 


VAa 2 + AS' 2t 




number 3. 


namc a 


Type b 


a(B1950) 


5(51950) 


[pc] 


q(B1950) 


<5(B1950) 


[PC] 


10 21 [cm- 2 ] 


[PC] 


comment 6 


224 


LMC N J0547-680 


I 


5h47.6m 


—68° 10' 


33 


5h48m 


-68° 10' 


55 


1.7 


79 




225 


LMC N J0547-704 


I 


5h47.8m 


—70° 42' 


87 


5h47.8m 


— 70°40' 




1.3 


28 


228 


226 


LMC N J0547-6953 


II 


5h48.7m 


—69° 54' 


61 


5h48m 


— 69°56' 


104 


1.6 


159 


223,227,270 


230 


LMC N J0555-681 


III 


5h55.9m 


—68° 10' 


73 


5h55.9m 


—68° 10' 




1.3 





271 


231 


T TV J /""I TV, T T n A A *~r />^rrt 

LMC N J0447-672 


I 


4h47.2m 


—67° 34' 


17 


4h47.ini 


—67° 36' 


190 


1.2 


28 


2,232 


232 


LMC N J0448-672 


I 


4h48m 


—67° 28' 


17 


4h47.2m 


—67° 18' 


82 


1.2 


229 


2 


233 


T TV If /"*! T\T T f\ A A t \ z^O"l 

LMC JN J 0449-681 


II 


4h49.8m 


—68° 16' 


17 


4h49.5m 


—68 12 


88 


1.3 


86 


6 


236 


LMC N J0458-661 


II 


4h58.7m 


-66°22' 


17 


4h58.4m 


-66° 24' 


75 


3.7 


83 


30,32,33,35 


237 


LMC N J0459-660 


II 


4h59.8m 


-66° 12' 


17 


5h0.2m 


-66° 12' 


* 


2.7 


70 




243 


LMC N J051 1-705 


I 


5hl2.3m 


-70° 54' 


17 


5hl2.7m 


-70° 54' 


8 


0.7 


87 




244 


LMC N J0511-692 


II 


5hl2m 


-69°32' 


17 


5hl2m 


-69°32' 


29 


1.1 





65 


247 


LMC N J0521-673 


I 


5h21.8m 


-67° 42' 


17 


5h22.ini 


-67° 36' 


34 


0.9 


117 




252 


LMC N J0528-672 


III 


5h28m 


-67°28' 


17 


5h28m 


-67° 28' 


29 


2.3 







253 


LMC N J0530-675 


I 


5h30.7m 


-68°0' 


17 


5h31m 


-67° 58' 


47 


2.1 


83 




258 


LMC N J0535-661 


II 


5h35.3m 


-66°20' 


17 


5h35.3m 


-66° 18' 


24 


0.7 


28 




266 


LMC N J0543-694 


II 


5h43.6m 


-69°48' 


17 


5h43.5m 


-69° 46' 


18 


3.6 


28 


207 


270 


LMC N J0547-695 


II 


5h48.1m 


-70°0' 


17 


5h48.6m 


-70°6' 


82 


2.0 


143 


226,227 


271 


LMC N J0553-682 


I 


5h53.3m 


-68° 24' 


17 


5h53.3m 


-68° 24' 


55 


1.9 








a Fukui et al. (2008) 
b Kawamura ct al. (2009) 

c position of peak integrated intensity. Data is from Fukui et al. (2008). 

d HI cloud radius defined as R = \J~^- Here S is the cloud area, calculated by summing the areas of all pixels detected above 80% of the peak integrated 
intensity level. Asterisks show HI clouds whose extent is poorly defined. 

C HI column density of position at peak integrated intensity estimated by using the relation; AT (HI) [cm -2 ] = 1.82 X 10 18 J* T(,(HI) dv [K km s -1 ]). 

'difference between CO peak position and HI peak position; (CO(a, <5)-HI(o, 5)) 

S HI clouds including two or more GMCs; The numbers show GMCs that are located in the same HI cloud. 




fthOQm 5100m 5h00m 4li30m 

Right Ascension (J2000) 

Fig. 1. — HI integrated intensity image (Kim et al. 2003) with contours of the CO integrated intensity 
(Fukui et al. 2008). The contour levels begin at 1.2 K km s _1 and have 3.6 K km s _1 intervals. 
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Fig. 2.— An example of the HI and CO line profiles at a(J2000) = 5 h 39 m 38 s , and <5(J2000) = -69°44'46". 
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Fig. 3. — The HI and CO datacubes that use for our analysis have somewhat different spatial and velocity 
resolutions; we have convolved both datasets to a spatial resolution of 40 pc x 40 pc, and a velocity resolution 
of x 1.7 km s _1 across the whole LMC. The pixels used in our analysis occupy the same volumes for both the 
HI and CO datasets. The 3cr detection limits for the CO and HI brightness are 0.21 K and 7.2 K respectively. 

(a) Histogram of the HI brightness temperature averaged within each 3-D pixel is shown in blue; the HI 
pixels with significant CO emission (i.e. > 3a) are marked in red. The width of the histogram bins is 10 K. 

(b) Ratio of the number of pixels with significant HI emission to those with significant CO emission within 
each 10 K bin. 




Fig. 4. — Histograms of the pixel-averaged HI brightness temperature where significant CO emission is 
detected for Type I (blue), Type II (yellow) and Type III (red) GMCs. Histograms are shown for the whole 
LMC, and for three different regions - Bar, North, and Arc - which are indicated in the right panel. 
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Fig. 5. — Velocity channel maps of the HI and CO emission associated with Type I, Type II, and Type 
III, GMCs. The grayscale images with black contours are HI brightness temperature (Kim et al. 2003) and 
the white contours are CO intensity (Fukui et al. 2008). The central velocity of the included channels is 
shown in the upper left of each panel. The lower right panel of each figure shows the total CO integrated 
intensity. The red contours and the red cross in the lower right panel of (a) and (b) present the 80 % level 
of the HI peak intensity and the position of the HI intensity peak, respectively, (a) GMC no. 225 in the 
NANTEN catalog. CO contour levels are from 0.105 K km s" 1 (3cr) with 0.105 K km s" 1 intervals. HI 
contour levels are from 20 K with 20 K intervals. This is a case where the 80 % level of the HI peak intensity 
is very extended relative to the size of the GMC. (b) GMC no. 134 and 135 between LMC4 and LMC5. CO 
integrated intensities are from 0.105 K km s _1 with 0.105 K km s^ 1 intervals. HI contour levels are from 20 
K with 20 K intervals. This is a case where the 80 % of the HI peak intensity provides a sensible definition 
of the HI envelope size, (c) GMC no. 192 and 202 with the HII region N159. CO integrated intensities are 
from 0.105 K km s _1 with 0.21 K km s _1 (6a) intervals. HI contour levels are from 20 K with 20 K intervals. 
This is a region where the HI emission is very complicated; it was not selected in the sample of 123 clouds 
that we use for our analysis. 
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Fig. 5(b). — continued 
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Fig. 5(c). — continued 
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Fig. 6. — Histograms of three parameters, (a) the distance from the kinematic center of the galaxy (kpc), 
(b) the CO linwidth (kim s _1 ) and (c) the CO cloud mass (M Q ) for the 123 GMCs used in our analysis. 
The solid lines show the distribution of these properties for all the clouds in the Fukui et al. (2008) catalog; 
the gray shaded regions represent the distribution for the 123 selected GMCs (see also Table 2). 
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Fig. 7. — The absolute difference between the central velocity of the HI and CO emission, |Vlsr(CO) — 
Vlsr(HI)|, versus the CO linewidth, (CO) for the 123 GMCs in our sample. The dotted line shows 
IVLsr(CO) - F LS r(HI)| - AV(CO). 




Fig. 8. — Histogram of the CO and HI linewidths toward the position of the CO peak position for the 
123 selected GMCs. The dotted lines show Gaussian fits to the histograms. The mean values (standard 
derivations) for the CO and HI distributions are 4.6 (1.6) and 14.1 (3.3) km s _1 respectively. 




Fig. 9. — Plot of the CO versus HI linewidth for the 123 GMCs in our sample. The red line is the regression 
line AV(CO) = (1.32 ± 0.04) + (0.23 ± 0.003)AV(HI), and the dotted line shows AV(HI) = 3 x AV(CO). 
The Spearman rank correlation coefficient is 0.39. 




VAa 2 + A5 2 [pc] 

Fig. 10. — Histogram of the projected separation between the CO and HI peak positions, as listed Table 
2. The separation between the CO and HI peak positions for four of the 123 clouds is greater than 300 pc; 
these clouds are not shown. 
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Fig. 11— The radius of the GMC, R(CO), versus the radius of the HI envelope, i?(HI) for the 62 
clouds with radii greater than 30 pc, as listed in Table 2. The red line is the regression line R(CO)= 
10~ a7±a5 i?(HI) 1 - 4±0 - 3 ; the Spearman rank correlation coefficient is 0.45. 
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Fig. 12. — The CO versus HI integrated intensity for the 123 clouds in our sample, averaged over the 
projected area of each GMC or HI envelope. Blue, green and red plots symbols represent Type I, Type II 
and Type III GMCs respectively. The black dotted line shows the regression line for all 123 clouds; /(CO) 
= 1CP 2 ' 6±0 2 /(HI) 11±01 . The regression lines and Spearman rank correlation coefficients (Ri, Rn, and Rm) 
for the three GMC types are shown at the top of the plot. 
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Fig. 13. — CO and HI emission in the N159 region. The grayscale image represents the combined 
ATCA+Parkes HI integrated intensity data (Kim et al. 2003). The black contours are integrated CO 
(J =3-2) emission by ASTE (Minamidani et al. 2008) and green contours are integrated CO (J =1-0) 
emission by NANTEN (Fukui et al. 2008). For CO (J =3-2), the contour levels are from 5 K km s" 1 in 10 
K km s _1 intervals; for CO (J =1-0), the contours levels are from 1.2 K km s _1 in steps of 2.4 K km s _1 . 
The thin black lines indicate the region observed in CO (J =3-2) emission by ASTE. The crosses indicate 
the positions of HII regions cataloged by Davies et al. (1976). 



